On Ising and dimer models in two and three dimensions 
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Motivated by recent interest in 2+1 dimensional quantum dimer models, we revisit Fisher's map- 
ping of two dimensional Ising models to hardcore dimer models. First, we note that the symmetry 
breaking transition of the ferromagetic Ising model maps onto a non-symmetry breaking transition 
in dimer language — instead it becomes a deconfinement transition for test monomers. Next, we 
introduce a modification of Fisher's mapping in which a second dimer model, also equivalent to 
the Ising model, is defined on a generically different lattice derived from the dual. In contrast to 
Fisher's original mapping, this enables us to reformulate frustrated Ising models as dimer models 
with positive weights and we illustrate this by providing a new solution of the fully frustrated Ising 
model on the square lattice. Finally, by means of the modified mapping we show that a large class of 
three-dimensional Ising models are precisely equivalent, in the time continuum limit, to particular 
quantum dimer models. As Ising models in three dimensions are dual to Ising gauge theories, this 
further yields an exact map between the latter and the quantum dimer models. The paramagnetic 
phase in Ising language maps onto a deconfined, topologically ordered phase in the dimer models. 
Using this set of ideas, we also construct an exactly soluble quantum eight vertex model. 
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I. INTRODUCTION 

Dimer models have long been of interest to statisti- 
cal mechanicians. ^^'^ In addition to their interest in vari- 
ous physical contexts, they have the striking feature that 
they are exactly soluble on any planar graph. ^ Following 
this insight. Fisher,^ after initial work by Stephenson,'' 
constructed a general mapping - reviewed below - from 
two-dimensional Ising onto dimer models, thereby relat- 
ing the solvability of the one to that of the other. In 
particular he related the partition function of the fer- 
romagnetic Ising model on the square lattice to that of 
the dimer model on the (now) Fisher lattice, which is 
sketched in Fig. 1. 

More recently, quantum dimer models (QDMs) have 
been formulated and studied.^ These live in Hilbert 
spaces spanned by dimer configurations of a given lat- 
tice and their Hamiltonians contain kinetic and potential 
energies that are naturally defined in this basis. These 
models were introduced to capture the dynamics of va- 
lence bond dominated phases of quantum antifcrromag- 
nets, with the particular intent of finding Anderson's hy- 
pothesized resonating valence bond (RVB) liquid^ - a 
hope realized recently.^" 

In this context, attention has been focussed on the na- 
ture of the quantum dimer phases, in particular their 
topological properties and low-energy excitations, and a 
recurring and useful theme has been their interpretation 
in gauge theoretic terms, particularly the identification 
of the RVB phase as a deconfined phase. ^^"^^ In Ref. 14 
this identification was made as an exact reduction of an 
"odd" Ising gauge theory to a quantum dimer model on 
the same lattice, in the extreme strong coupling limit. 
As Ising gauge theories in d = 2 -I- 1 are dual to Ising 
models, this allows dimer models to be exactly related 
to frustrated quantum (transverse field) Ising models in 
d = 2 -f 1 in the limit of weak transverse fields and 
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FIG. 1: Left: the Fisher lattice, a decorated square lattice. 
The six sites in the basis correspond to one square lattice 
site. A dimer model on this lattice can be used to calculate 
the partition function of the ferromagnetic Ising model on the 
(direct) square lattice. Right: The frustrated Fisher lattice, 
suitable for a dimer model representing the Ising model on 
the fully frustrated (dual) square lattice. 



thence to a set of ideas and techniques for obtaining their 
phase diagrams. Very recently, Misguich et al. have con- 
structed an exactly soluble dimer model on the kagome 
lattice. This model maps onto a transverse field Ising 
model with zero exchange, thereby allowing the entire 
spectrum and correlations to be determined. 

In this note, we further explore the above connections 
between Ising and dimer models in two and three dimen- 
sions. We begin by reviewing Fisher's construction in 
two dimensions which utilises the loop model generated 
by the Ising high-temperature expansion, which in turn 
is mapped onto a dimer model on a decorated lattice. 
We note that the symmetry breaking transition in the 
two dimensional ferromagnetic Ising model maps onto a 
non-symmetry breaking, deconfinement transition of test 
monomers in the dimer model. As an aside we point out 
that the dimer formulation provides an immediate insight 
into how the Ising model yields a lattice theory without 
doubled fermions. We next observe that in the presence 
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FIG. 2: The mappings of the Ising model onto dimer models. 



of frustration in the Ising model, Fisher's construction 
leads to negative weights in the dimer model. To remedy 
this, we introduce a modification of Fisher's construction 
which proceeds via an intermediate map to a generalised 
(non-hardcore) dimer model on the dual lattice which is 
then decorated to produce the hardcore constraint. As 
an example, we solve the fully frustrated Ising model on 
the square lattice as a dimer model on a modified Fisher 
lattice. 

Finally, we turn to a class of three-dimensional clas- 
sical Ising models — these are general (frustrated or un- 
frustrated) two-dimensional Ising magnets stacked ferro- 
or antiferromagnetically in the third direction, which can 
therefore be mapped onto transverse field Ising magnets 
in 2-I-I dimensions by taking a continuum limit in the 
third direction. Armed with the modified Fisher con- 
struction we are able to map these onto quantum dimer 
models where the quantum dynamics induced by the 
transverse field translates to a particular "resonance" dy- 
namics of the kind studied by Misguich et al. This is 
particularly interesting from the viewpoint of the dimer 
models since it allows known lore on Ising models to 
be transcribed into statements about the former. Quite 
generally, the paramagnetic phases of the Ising models 
map to deconfined, topologically ordered phases in the 
dimer models. Specifically, for the case of the ferromag- 
netic Ising model we find a dimer transition character- 
ized solely by confinement and the loss of topological or- 
der. We also observe that the dimer models are exactly 
equivalent to Ising gauge theories, as they must be since 
the latter are dual to Ising models in 2-|-l dimensions. 
Consequently their spectrum in the deconfined phase is 
understood in terms of Ising vortices or visons. The rela- 
tionships between these mappings are displayed in Fig. 2. 

We close with some remarks concerning possible exten- 
sions of this work, including the construction of an ex- 
actly soluble quantum eight- vertex model along the lines 
of Ref. 15. Such a model may also be of interest in the 
context of orbital current models, such as the d-density 
wave models considered in Ref. 17. 



II. THE ISING TRANSITION IN d = 2 AS A 
DECONFINEMENT TRANSITION 

In this and further sections, we consider classical Ising 
models defined by the Hamiltonian 



iJ = - ^ Jij(Ji<Tj, 



(2.1) 



where cr^ is a Pauli matrix, and Jij is the strength of bond 
(ij). The sites of the lattice, labelled by « = 1---N, 
together with the bonds (ij) with nonzero Jij form a 
graph, Q, the interaction graph. In the following, we will 
use the following definitions: /3 = l/fc^T; Kij = [3 Jij] 
Vij = 1/ Wij = tanh Kij . Although we have formally writ- 
ten a quantum Hamiltonian, the problem is, of course, 
still classical. 

We first briefly review Fisher's mapping.^ Its starting 
point is the high temperature expansion of the partition 
function, ZTrcxp(— /3iJ), of the Ising model: 



Z = 2 



N 



Y\ cxp{-Kij) cosh{Kij) 



T(%;t/). (2.2) 



Here, T(vij; Q) is the crucial quantity: it is the sum over 
all loop coverings, labelled by T{Q), of the graph Q. These 
loops can intersect one another, provided an even num- 
ber of links emanate from each site. The weighting of a 
particular covering is given by the product of Vij of all 
the links of its loops, so that 



'^{vij;g) = E n 



(2.3) 



Fisher's mapping turns the resulting loop model on the 
graph g into a dimer model on a decorated graph, V, by a 
decoration procedure outlined in Ref. 6, see in particular 
its Fig. 6 for the general decoration rule and Fig. 7 for 
the explicit example of the square lattice magnet. The 
resulting lattice is depcited in Fig. 1. There are now two 
types of bonds, the original ('external') ones and the new 
('interior') ones. 

The crucial property of Fisher's mapping is that, if 
g is a planar graph, so is V. The partition function 
of the dimer model on V, Z-p, which is related to that 
of the original Ising model in a simple way, can in that 
case be evaluated by Pfaflian methods; by Kasteleyn's 
theorem,^ one can assign directions to each bond of V 
such that the product of orientations when traversing 
any elementary plaquette in a clockwise direction is odd. 
These orientations can be used to define a matrix. A, 
which has entries Aij — -|-1(— I) if bond (ij) is oriented 
for i to j (from j to i). With fugacities of the external 
and internal bonds given by Wij and 1, respectively, one 
has Zt, = Pf A. 

A nonanalyticity in Z bequeathes one to Z-p. For in- 
stance, the Fisher lattice dimer model obtained from the 



square lattice ferromagnetic Ising model must exhibit a 
phase transition as the ratio of external to internal fu- 
gaoities is varied. Nevertheless the Ising character of the 
transition seems to have disappeared en route — there no 
longer is an Ising symmetry to break. In the follow- 
ing paragraphs, we will see that the dimer model does 
not break any symmetries whatsoever. Instead the Ising 
transition has turned into a confinement-deconfinement 
transition in which the free energy to separate two test 
monomers is the appropriate diagnostic. The universal- 
ity classes of the two transitions coincide, as they must, 
by virtue of both of them mapping onto the theory of a 
single gapless Majorana Fermion. 

As V is periodic, one can find the partition func- 
tion by Fourier transformation; Z = Pi A = v det A = 

y^rig det A{q), where A{q) is the Fourier transform of A 
at wavevector q. The specific free energy, !F, is given by 

-AT/jjF = lnZv = (l/2)^ln(deti(g)) . (2.4) 

This expression does indeed reproduce Onsager's for- 
mula for the square lattice Ising model for any cou- 
pling strength K.^'^^ In particular, nonanalyticities occur 
when det A{q) = 0. With^ 

det A{q) = (1 -t- w^)^ — 2w{'w^ — l)(cosQ'a; + cosg^), 

this happens only for g^. = = and w = Wc = 1 + \/2 
so that Kc = ln(l -|- •\/2)/2, as it should. 

What is the nature of the phases on either side of 
u>c? For w oo, only one dimer configuration survives, 
namely one in which all external bonds and the bond 
linking the two internal sites are occupied. This configu- 
ration breaks no lattice symmetries, and is in that sense 
not a crystal. In the opposite limit, w = 1, all dimer 
configurations have equal weight and their ensemble also 
respects all lattice symmetries. Evidently the transition 
does not involve symmetry breaking. 

Instead, the two phases differ in their response to the 
insertion of a pair of test monomers (sites not part of 
a dimer). It is not hard to see that the high tempera- 
ture phase is confining. If one places a monomer on one 
end of an external leg, the site on its other end has to 
pair up with another site in its cluster (group of sites 
obtained from one original site by decoration), which in 
turn leaves the partner of the latter site unpaired. Two 
monomers placed a distance L apart therefore exact a 
cost in free energy proportional to the minimal number 
of unoccupied external bonds, which is proportional to 
L. This point is thus in a confined phase. Intuitively, 
the low temperature phase involves a true dimer fluid 
which should therefore allow monomers to be separted 
with finite free energy cost. 

These statements can be made precise by tracking the 
spin-spin correlation, {(JiCTj), from the spin formulation 
into the dimer formulation. The spin correlator can, in 
the context of the high-temperature expansion, be ex- 
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FIG. 3: The open string configurations (four rightmost plots) 
are generated with a monomer (filled circle) placed on the up- 
per interior point of the cluster. The other monomer location 
(leftmost plot) generates the remaining four string configura- 
tions, obtained by a reflection of the four displayed ones along 
a diagonal. 



pressed as a loop model on the square lattice contain- 
ing, along with the closed loops, one open string run- 
ning between the sites i and j.^^ This loop-string model 
can in turn be cast in terms of a dimer model with 
two monomers (Fig. 3). The vertices at the ends of the 
string can be encoded in the monomer-dimer model by 
summing over four partition functions, Zmd, with the 
monomers in the clusters i,j being placed independently 
onto one of two interior points of the cluster. The pari- 
tion function of the loop-string model is thus given by the 
sum over the four monomer-dimer partition functions: 
Zis{i,j) = Y.Zmd{i,3)- 

Next, one uses {uiOj) — Zis{i,3)IZD- It follows that 
in the high temperature (paramagnetic) phase, where 
{oiUj) vanishes exponentially, all the monomer pairs are 
confined. Conversely, in the low temperature phase, 
{(JiCTj) decays to a constant, which implies that at least 
one (and probably all) monomer pairs are deconfined. 

Two comments on related issues are worth making. 
First, dimer models generally exhibit topologically dis- 
connected sectors under local moves a classical pre- 
cursor of the notion of topological order for quantum 
problems^^ which we shall invoke below for quantum 
dimer models. On the Fisher lattice there are four such 
sectors on a torus corresponding to the various combi- 
nations of an even or odd number of dimers intersecting 
the two non-trivial loops. The Fisher map represents the 
Ising partition function as a sum over all four sectors 
of the dimer model. Evidently, all Ising spin configura- 
tions in a finite system are accessible from one another 
by local spin fiips, and the topological distinction in the 
dimer model is an artefact of the bookkeeping from the 
perspective of the Ising model. 

Second, dimer problems are solved by lattice Majorana 
fermions. As the unit cell of the Fisher lattice has 6 sites, 
microscopically the dimer formulation leads to as many 
Majoranas. However, at the transition fugacity Wc, two 
eigenvalues of A{Q) vanish. These two combine to form 
one single Majorana Fermion, to yield the known crit- 
ical theory as in the Ising model. We note that given 
the lore on fermion doubling,^" it is somewhat surpris- 
ing that we obtain only a single Majorana Fermion from 
our lattice problem. While this has been commented on 
from a different persepctive before, we note that in the 
dimer formulation this conclusion arises from the mani- 
festly asymmetric nature of our lattice, no longer invari- 
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FIG. 4: Decoration of vertices of coordination q (top) into 
clusters to generate the lattice on which a frustrated (mid- 
dle) or unfrustrated (bottom) Ising model is represented by a 
hardcore dimer model with positive fugacities. 

ant under reflections along both x and y axes. This leads 
to lattice derivatives different from the ones normally en- 
countered when discretising the continuum. 

III. FRUSTRATED ISING MODELS AS DIMER 
MODELS WITH POSITIVE FUGACITIES 

A. The modified Fisher construction 

Fisher's construction is quite general. The high tem- 
perature expansion generates a loop model on any inter- 
action graph and the rules for generating the clusters 
out of the vertices depend only on the number of legs of 
each given vertex. 

One prehaps not so desirable feature of the Fisher 
mapping is that it leads to negative dimer fugacities 
in the case of frustrated models: there, not all inter- 
actions can be chosen to be ferromagnetic, and hence 
the Vij = 1/wij = tanh Kij are negative for some of the 
bonds. 

Here, we present a way of addressing this problem via 
a modified Fisher construction on the dual of a two- 
dimensional planar lattice. It proceeds by mapping each 
Ising spin configuration onto a link configuration, {r}, 
on the dual lattice. 

A given link of the dual lattice is occupied (t = 1) if 
and only if the bond of the direct lattice it crosses is frus- 
trated. Such a link configuration has the property that 
the site of the dual lattice at the centre of a frustrated 
plaquette has an odd number of occupied links emanat- 
ing from it; for an unfrustrated plaquette, this quantity is 
even. Clearly a spin configuration and its Ising reversed 
counterpart map onto the same link configuration on the 
dual lattice. 

Such soft dimer configurations on the dual lattice can 

be converted into hardcore dimer coverings by suitably 
decorating each site. Depending on the even or oddness 



of the number of occupied links at the site, two different 
decoration operations arc required, leading to two dif- 
ferent types of clusters. The decoration transformation 
for odd sites is shown in Fig. 4. For even sites, one re- 
quires a slightly different prescription as in Fisher's orig- 
inal construction, as that included interchanging empty 
and occupied links. The corresponding construction is 
also given in Fig. 4. 

The dimer fugacities arc determined by the Boltzmann 
factors which come with the presence of a frustrated 
bond. We find it convenient to add a constant term to 
the Hamiltonian so that an unfrustrated bond has energy 
and a frustrated bond has energy 2 J, so that the dimer 
fugacity link (ij) is given by w = exp{—2Kij), which is 
always positive and between and 1. As in Fisher's orig- 
inal construction, Ising models on planar lattices lead to 
planar dimer models, which can hence be solved using 
Kasteleyn's theorem. 

Three comments on the differences between the Fisher 
construction and the modified Fisher construction are 
in order. First, as is appropriate for a dual construc- 
tion, high and low temperatures trade places when we 
compare the original Fisher construction to the modified 
construction — e.g. the equal fugacity dimer model corre- 
sponds to zero and infinite temperature respectively. Sec- 
ond, in the modified construction the spin model maps 
onto a single topological sector of the dimer model. The 
remaining sectors are generated by considering different 
boundary conditions for the spin model. Third, the orig- 
inal mapping does not relate individual spin configura- 
tions to dimer configurations but the modified mapping 
does, upto a twofold ambiguity coming from global Ising 
reversal. This will be important in making a connection 
between Ising models and quantum dimer models below. 
But first we apply the modified construction to the solu- 
tion of a classical frustrated Ising model. 



B. The fully-frustrated Ising model as a dimer 
model 

As an an illustration of the above technique, we now 
use it to solve Villain's odd model, also known as the 
fully frustrated Ising model on the square lattice. This 
model is defined for spins on the square lattice with near- 
est neighbour interactions such that each plaquette has 
an odd number of antiferromagnetic interactions. This 
can for example be achieved by choosing all horizontal 
bonds to be ferromagnetic and the vertical bonds in fer- 
romagnetic rows alternating with antiferromagnetic ones. 
Transcribed to the dual lattice, which is again square, 
this corresponds to an odd number of occupied links em- 
anating out of each plaquette. The corresponding deco- 
rated lattice (Fig. 1) has five sites per cluster. 

The Kasteleyn matrix on this lattice can be partially 
diagonalized by means of a Fourier transform. While 
the lattice itself has five sites per unit cell (see Fig. 5), 
choosing appropriate signs for the matrix requires that 
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we consider a doubled unit cell in one direction which 
we take to be the horizontal direction. This yields the 
10 X 10 Fourier transformed matrix, A{q), 



FIG. 5: Arrow convention for the Kasteleyn matrix of the 
frustrated Fisher lattice. 



A{q) = 
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Its determinant, detA{q), is given by 

deti(g) =2^2 (2(1 + + (-1 + ^2)^ 



(3.1) 



cos(g'x) - (-1 + u^)^ cos(2g'y)j , 



(3.2) 



which yields the dimensionless free energy per site, 

The first term is the ground state energy, as there is one 
frustrated bond for each pair of sites. At zero temper- 
ature, the second term gives the groimd state entropy, 
which integrates to the well-known result G/tt, where G 
is Catalan's constant. ^'^'^ After a little algebra, our ex- 
pression Eq. 3.3 can be shown to agree with Villain's 
result for the partition function given in Appendix 2 of 
Ref. 22 at all temperatures. 

A direct examination of Eq. 3.3 shows that the model is 
critical only at T = 0, where u = 0. As a final exccrcisc 
we will now compute the divergence of the correlation 
length as T — > by means of our solution. 

To find its behaviour in the critical regime, we con- 
sider the case of small u. Near u = 0, we can expand 
det(^(g'))/(2u^) to find that it varies as 

2 + cos(gx) — cos(2gj^) + 2'U^(2 — cos{qx) + cos{2qy)X3A) 
In particular, at it = 0, it vanishes at pf'' = (tt, 0) and 

To find the correlation length, it is necessary to com- 
pute the Green function by inverting A. This is done 



by inverting A{q) to obtain and then carrying out the 
inverse Fourier transform on G{q) = A~^{q). If one is 
only interested in the correlation length and not the de- 
tails of the correlations, it suffices to do the Fourier inte- 
gral jj cPqexp(iqr)G{q) asymptotically, using the prop- 
erty that due to the inversion process, the structure of 
G{q) = g{q)/ det A{q), where g{q) denotes a cofactor. We 
thus have to do integrals of the type 

G - JJ d^g exp(ig • fR)/f{q) , (3.5) 

where the zeroes of det A{q) determine the locations of 

the poles of the integrand - we have checked that the 
cofactors do not all vanish at the locations of the poles. In 
this equation, we have emphasised the two-dimensional 
nature oi q = {qx,qy) by writing it as q for the time 
being, f = (cos 9, sin 9) is a unit vector in the direction 
of which the correlations arc to be computed, and R 
00 is the quantity in which we will evaluate the integral 
asymptotically, so that {x,y) = Rr. 

Let us first carry out the integral over q^. 

Ix = j dqxex.p{Riq^ cos 9)/ f{q). (3.6) 

From Eqs. 3.2, 3.4, we find that the location of the poles 
can be written in terms of g* = {2q'^,q'y) so that pi = 
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ffl'^ + (f , where we have inserted a factor of 2 in front of 
q'^ to make the square symmetry more apparent: 



?'2 



/2 , ,2 



= -Au^ = -a 



(3.7) 



so that the poles he at q'^ = ±^ a^v? + q'y^. Depend- 
ing on the sign of cos 6, we therefore choose an integra- 
tion contour that runs through the following points in 
the complex q^ plane: (0, 0) (27r, 0) (2?:, ±oo) 
(0, ±oo) — > (0,0). Due to the periodicity of the inte- 
grand, the contribution from the vertical contours cancel. 
The contribution from the horizontal contour at ±7c» 
vanishes for i? ^ oo as exp(— i?| cos 6*1), so that only 
picks up a contribution from the enclosed poles, where 
we denote the residue as 1//'^ [Pi (?(,)] • 
To find G, we then need to determine 



G 



/■2tt 

/ dqyl^ 
Jo 



(3i 



X exp 



rlPliq'y)] r[P2{q'y)] 



-Ri \ cos 0\ ^ a^u^ + q'y^ — iq'y sin ( 



This integral can be treated asymptotically using 
Laplace's method. The dominant contribution comes 
from points where the argument of the exponential is 
stationary: 



iq'sm.e\=Q. (3.9) 



One finds q'y = iausva.6, so that at large distances, G de- 
cays exponentially as exp(— i?/^), where the correlation 
length = 2u. As T ^ 0, the correlation length of the 
Ising model therefore displays a divergence proportional 
to exp(2J/A;BT), in agreement with the result given in 
Ref. 23, which was obtained using the transfer matrix 
technique. 



IV. THE THREE DIMENSIONAL ISING 
MODEL AS A TWO DIMENSIONAL QUANTUM 
DIMER MODEL 

In this section we use the modified Fisher construction 
to reformulate a class of three dimensional Ising mod- 
els as quantum dimer models in 2 -f- 1 dimensions. The 
interest of this mapping is that it yields non-trivial in- 
formation on the quantum dimer models by transcribing 
known lore on the Ising models — no advances in solubil- 
ity are entailed. 

The prototypical example of the Ising models of inter- 
est is the nearest neighbour ferromagnetic Ising model 
on the cubic lattice. As is well known, this is equiva- 
lent to a 2 -|- 1 dimensional transverse field Ising model 
in the "time continuum" limit. Briefly, this proceeds by 




FIG. 6: One kinetic energy term of the quantum dimer Hamil- 

tonian generated by the simple cubic Ising ferromagnct has a 
matrix element between the two configuration shown here. 



the recognition that the phase structure of the model 

is unchanged if we take the cubic lattice to be a set 
of stacked square lattices and allow unequal couplings 
in the plane and in the stacking direction, H{Jt,Js) = 

-JrY^u,x(Tfcr^ - JsY.Ui)^i^^j where the sums {ij) 



and {ij) run over nearest neighbour sites in the stack- 
ing and planar directions. The anisotropic scaling limit 
exp(2_ft"^) GO can be identified with the Trotter-Suzuki 
decomposition of the imaginary time path integral of 
the two-dimensional transverse field Ising model with the 
quantum Hamiltonian H — ~JsJ2{ij) '^i'^j ~ ■ 

The Hilbert space of this quantum model is spanned by 
all classical Ising configurations of the two dimensional 
classical model and the transverse field moves the system 
between these configurations. As the classical configura- 
tions can be related to dimer configurations on the Fisher 
lattice by the modified Fisher construction, it should be 
intuitively clear that the model can equally well be cast 
as a resonance dynamics in the space of dimer configura- 
tions. 

The quantum dynamics induced in this way can be vi- 
sualised by noting that the effect of the transverse field is 
to fiip individual spins. In the modified Fisher construc- 
tion, this corresponds to replacing empty external links 
surrounding the spin by occupied ones, and vice versa. In 
the dimer model, this also entails moving internal bonds 
to accomodate the change in external bonds. Here, it 
is important to note that the configuration of internal 
bonds of a cluster, given a set of external bonds, is unique 
and determined only by the by the external bonds belong- 
ing to that cluster. The quantum dimer Hamiltonian is 
therefore strictly local, although it includes a number of 
kinetic energy terms and loop fiips of varying length with 
exactly the same strength — all nonzero off-diagonal ma- 
trix elements equal — F. An example of a dimer move 
present in the current Hamiltonian is given in Fig. 6. 
In addition the Ising nearest neighbour interaction — the 
number of frustrated bonds translates into a potential 
energy for dimers on the external bonds. 

The resulting quantum dimer Hamiltonian of this 
model can be schematically written as 

Bu = V + f 

= -E^^i^5+^E^^^ (4-1) 

(ij) n □ 

The Jij of internal bonds are zero and the sum 
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runs over all closed loops made up of the external bonds 
surrounding a spin of the direct lattice and any of the 
bonds of the clusters on which these bonds terminate. 

stands for alternating raising and lowering opera- 
tors, • • • T~f.T^iTjl^ • • • as one goes around the loop; this 
form preserves the hardcore condition. 

In the above we have ignored one subtlety, namely that 
the map between Ising configurations and dimer configu- 
rations is two-to-one and the Ising dynamics does connect 
a state and its Ising reversed counterpart. The solution 
to this lies in considering combinations of a given state 
and its reversed counterpart^^ 

|{a0)e/o = ^(|{aj)±|{-aj)) (4.2) 

that are even and odd under global Ising reversal — a 
property respected by the Ising dynamics. In order to 
fix the "up" states we choose those so that a particular 
spin is always up in them. Now both sets of states can 
be mapped onto dimer states but the dimer dynamics is 
slightly (but importantly) different in the two sectors. In 
the even sector it is exactly what we described above. 
In the odd sector most matrix elements are the same 
but the ones that involve the chosen spin acquire an ex- 
tra minus sign. Since F can be chosen negative without 
loss of generality, by the Pcrron-Frobenius theorem the 
ground state will always be in the even sector, although 
the first excited state need not be. So for the purposes of 
determing phase structure one can ignore this complica- 
tion entirely. From the perspective of the dimer model, 
which is what we will take in the remaining, its dynamics 
will be represented solely by the even states of the Ising 
model. 

At zero temperature, there is one parameter in this 

problem, namely the ratio of transverse field to bond 
strength, F/J^ and two phases that meet at a critical 
point. In Ising language, the two phases are, of course, 
the ferromagnet and paramagnet. As the reader may an- 
ticipate from our previous classical considerations, the 
dimer transition is between a deconfining phase at large; 
F/ Jg and a confining phase at small T/Jg neither of which 
break any lattice symmetries. 

The dcconfined phase is of particular interest. Its 
point F/Jg = oo is the "Rokhsar-Kivelson" point of the 
model where the ground state wavefunction is the equal 
amplitude sum over all dimer configurations, which is 
the canonical short ranged RVB state. As discussed in 
Ref. 15, the elementary Ising excitation is a a spin an- 
tialigned with the transverse field with energy 2F, and 
the entire spectrum is composed of them. In the dimer 
model the fiipped spin translates into a vortex in which 
dimer configurations pick up a minus sign if the number 
of dimers on a string extending from the chosen plaquctte 
to infinity is odd (labelled by dotted bonds in Fig 7). The 
dimer model is exactly equivalent to the standard Ising 
gauge theory, and the vortex is therefore an Ising vortex 
(vison). For a system on a torus, these vortices need to 
be created in pairs, so that the minimal excitation energy 




FIG. 7: A spin antialigned with the transverse field at the 

centre of the plaquottc denoted by the cross translates into a 
vortex excitation in the quantum dimer model. 

is in fact AV}^ 

The Ising states with a single flipped spin have disap- 
peared in the course of the duality transformation. In 
fact, together with the other states with an odd num- 
ber of fiipped spins, they make up the sector, odd under 
global Ising reversal, which was discarded en route to the 
Ising model. The ground state at the Rokhsar-Kivelson 
point in this sector is the first excited state of the Ising 
model and has a degeneracy proportional to the system 
size. The equal amplitude superposition of the configu- 
rations in that sector corresponds to the reference spin 
being antialigned with the field. 

As observed before there are four dimer sectors on the 
torus (and 4^ on lattices of genus g) which cannot be 
connected by local dimer moves. As these sectors corre- 
spond to different boundary conditions in the Ising prob- 
lem, and the latter have only an 0{e~^^^) effect on the 
ground state energy in the paramagnetic phase (^ is the 
correlation length), it follows that they lead to four expo- 
nentially degenerate ground states. Hence the dcconfined 
phase displays topological order in the sense of Wcn.^^ 

Finally we observe that the dimer model on the Fisher 
lattice is exactly equivalent to the standard Ising gauge 
theoy on the square lattice - dubbed the "even" Ising 
gauge theory in Ref. 14 due to the nature of its con- 
straint. That constraint states, in the Hamiltonian for- 
mulation, that the number of units of Ising electric flux 
entering/leaving a site (the distinction is irrelevant for 
Ising variables) must be even. If we identify a dimer on 
the external bonds of the Fisher lattice with such a flux, 
we recover this constraint. This conclusion is only to be 
expected since the even Ising gauge theory is dual to the 
transverse field Ising model. 

Interestingly, the dimer constraint on the Fisher lat- 
tice, that the number of dimers coming out of a given 
site be one, is a U{1) constraint. As discussed in Ref. 14, 
this gives the dimer model the character of a C/(l) gauge 
theory at the lattice scale. As in the case of the trian- 
gular lattice QDM, the topological sectors indicate a low 
energy structure with Ising character - which would be 
the general expectation for a dcconfined (RVB) phase. 
What is special here is that, unlike the triangular QDM, 
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wc arc able to put all states of the Fisher lattice QDM 
under consideration into correspondence with states in 
the Ising gauge theory. 

Our basic considerations in this section can be gener- 
alized to a wide class of Ising models that permit a time 
continuum limit to be taken. These include models that 
can be viewed as identical planes stacked ferromagneti- 
cally or antiferromagnetically (or in any alternating com- 
bination). Depending on the frustration in the planes, we 
will obtain a QDM on a specific lattice obtained by the 
decoration procedure. For example, the FFIM stacked 
ferromagnctically will give rise to a QDM on the lattice 
shown in Fig. 1, which will be exactly equal to the "odd" 
Ising gauge theory in which the number of Ising fluxes 
leaving a site is odd. This QDM will exhibit a transi- 
tion from a deconfined phase to a confining phase that 
does break lattice symmetries, as first discussed in the 
language of Ising models in Ref. 27. 



V. COMMENTS ON THE SOLUBLE KAGOME 
QUANTUM DIMER MODEL 

The dimcr model we obtain at F/ J, oo has much 
in common with the kagome dimer model discussed by 
Misguich et al.: their quantum dimer model also has a 
range of kinetic terms around a given hexagonal plaque- 
tte, and it displays an analogous excitation spectrum con- 
sisting of Ising vortices. Nevertheless, the solubility and 
beautiful simplicity of their model does not derive from a 
Fisher construction. Instead, their problem has another 
very useful ingredient, namely the mapping of hardcore 
dimers to arrows on the kagome lattice, and indeed any 
other lattice consisting of corner-sharing triangles. 

Not until they have executed this mapping docs our 
discussion of the Rokhsar-Kivelson point parallel their 
model, as one can consider the arrows as bond variables of 
a triangular lattice Ising model (or cquivalcntly, link vari- 
ables of a honeycomb gauge theory) which, crucially, 
has zero exchange strength. The quantum dynamics they 
study consists of the Wilson loop operation Y\q t^. 



A. An exactly soluble quantum eight vertex model 

To illustrate this correspondence, and because it is of 
interest to ask if the arrow mapping used by Misguich et 
al. can be generalised, we note that their approach can 
be used to define an exactly soluble quantum eight vertex 
model. This model can also be given an interpretation in 
terms of singlet bonds. 

The basic geometric object of this model are squares 
(instead of triangles), which arc again arranged to share 
corners. At the corners of the squares reside Ising degrees 
of freedom (arrows), which point either in or out. In 
addition, wc impose the constraint that an even number 
of arrows point out. 
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FIG. 8: A lattice of corner sharing squares (thick lines), which 
is the medial lattice of the lattice denoted by thin lines, which 
is in turn dual to the lattice denoted by dashed lines. 

The simplest case occurs when the squares (denoted 
by fat lines in Fig. 8) are arranged on a square lattice 
(thin lines) so that they share corners. The different 
arrow configurations can be mapped onto configurations 
of the eight vertex model on the dual square lattice by 
identifying the arrows (which by construction live on the 
links of the dual lattice) with the arrows of the vertices 
of the eight- vertex model. 

One can further interpret the arrows as link variables 
of an Ising model. The sites defining this model live on a 
further dual lattice (dashed lines). For the case of a bi- 
partite thin lattice, one can label a bond frustrated if the 
arrow points from sublattice A to sublatticc B and un- 
frustrated otherwise. (In the case of a non-bipartite thin 
lattice, one can choose a reference arrow and a reference 
spin configuration and identify the two). The constraint 
of the eight-vertex model then becomes a constraint on 
the product of exchanges, Hn Jij = 1; around a plaquette 
of the dashed lattice. However, this only imposes a con- 
straint on allowed states - the strength of the exchanges 
vanishes. From this method, one can easily visualise the 
known result that the number of eight-vertex configu- 
rations on the square lattice equals 2^, as there are no 
constraints on the Ising model on the dashed lattice. This 
result is straightforwardly generalised to any lattice de- 
fined by the midpoints of corner-sharing squares, as ex- 
plained in Ref. 15 for lattices defined by midpoints of 
corner-sharing triangles. 

This eight-vertex model can again be endowed with a 
quantum dynamics generated by a Wilson loop action on 
the links of the elementary plaquette of the thin square 
lattice. This becomes the zero-exchange transverse field 
Ising model on the dashed square lattice. The results, 
such as on the gap, 4r, or on the ultra-short correlation, 
follow just as they did before.^''' 

A possible interpretation of the different eight-vertex 
configurations in terms of singlet bonds is as follows (see 
Fig. 9). Let spins = 1/2 reside on the vertices of the 
fat square lattice. If no arrows point into a given square. 
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FIG. 9: Possible interpretation of the eight vertex configu- 
rations (denoted by arrows on thin lines) in terms of singlet 
bonds (denoted by dimers) between spins 1/2 on the corners 
of the square. The rightmost figure denotes an overall singlet 
of the four spins, of which there are two linearly independent 
ones. 

there are no singlet bonds between any pair of spins of the 
square. If two point in, there is a singlet bond between 
the two. If four point in, the spins on the square form 
some collective singlet state. The choice of equal fugaci- 
ties for all vertices amounts to disregarding the entropic 
contribution of the two linearly independent collective 
singlets. 

As a related (and known) combinatorial result, we no- 
tice in passing that the number of non-overlapping (but 
possibly intersecting) loop configuration on a lattice can 
be trivially evaluated in the same spirit. One can think of 
the loops as domain walls of an Ising model on the dual 
lattice. Up to a factor 2 from the global Ising symmetry, 
the number of domain wall configurations is then given 
by the size of the configuration space of the Ising model, 



which is of size 2 for a dual lattice with N sites. 



VI. SUMMARY 

Two dimensional Ising models are intimately related 
to dimer models. By means of the modified Fisher con- 
struction introduced in this paper, this connection can 
be exhibited configuration by configuration. The Ising 
transition in two dimensions maps onto a deconfinement 
transition in the dimer language as does the Ising transi- 
tion in three dimensions which now appears in a quantum 
dimer model. This last equivalence provides an instruc- 
tive example of a topologically ordered RVB phase and a 
transition from it to a confining phase, and one for which 
all the details are known. These ideas can be extended 
straightforwardly to generate an exactly soluble quantum 
eight vertex model. 
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